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We consider the phase stability of a local oscillator (or laser) locked to a cavity QED system 
comprised of atoms with an ultra-narrow optical transition. The atoms are cooled to millikelvin 
temperatures and then released into the optical cavity. Although the atomic motion introduces 
Doppler broadening, the standing wave nature of the cavity causes saturated absorption features to 
appear, which are much narrower than the Doppler width. These features can be used to achieve 
an extremely high degree of phase stabilization, competitive with the current state-of-the-art. Fur¬ 
thermore, the inhomogeneity introduced by finite atomic velocities can cause optical bistability to 
disappear, resulting in no regions of dynamic instability and thus enabling a new regime accessible 
to experiments where optimum stabilization may be achieved. 


INTRODUCTION 


Today’s ultra-precise and accurate atomic clocks con¬ 
tinue to make important contributions to fundamental 
physics as well as applied technology. Atomic clocks have 
imposed significant constraints on the drift of fundamen¬ 
tal constants [Il-Q, may have the potential to enhance 
the sensitivity of gravitational wave detectors, and have 
provided ultimate tests of the general theory of relativity 
[^. With the current stability of optical clocks at the 
1x10“^® level, there are prospects for applying atomic 
clocks for the detailed mapping of the Earth’s gravity 
field Hi. 

A highly stabilized laser is an integral component 
of high precision measurements, such as optical atomic 
clocks and precision spectroscopy. Current technology 
for achieving highly phase stable laser sources relies on 
locking a laser to a high-Q reference ULE glass cavity 0- 
Q. The phase stability of this method is currently limited 
by the thermal noise induced in the mirrors, spacers, and 
coatings of the reference cavity , but has been signif¬ 
icantly reduced over the past few years with new engi¬ 
neered materials [Hllll- 

As an alternative approach to overcoming the thermal 
noise problem, it was recently proposed [ISj to lock the 
laser to the saturated resonance feature exhibited by a 
collection of atoms with an ultra-narrow electronic tran¬ 
sition trapped in an optical cavity. Here, the atoms were 
assumed to be trapped in an optical lattice inside the 
cavity. Due to the narrow atomic line, such a system 
would typically operate in the parameter region corre¬ 
sponding to the bad cavity limit of cavity QED. There, 
the atomic linewidth is significantly narrower than the 
cavity linewidth. In contrast to the reference cavity stabi¬ 
lization method described above, the cavity QED method 
offers a distinct advantage since no drift compensation is 
needed. 

The cavity QED system exhibits optical bistability in 
the intracavity intensity Bui where several solutions 
exist for the steady-state intracavity field. Working at an 
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FIG. 1. (color online). Schematic of the cavity QED experi¬ 
ment with a thermal sample of atoms with Doppler width Fd, 
each of which have a narrow optical transition of width 7a << 
Fd. The atoms are probed with a carrier with frequency oJc, 
and two sidebands located at tJc±FSR, where FSR is the free 
spectral range of the cavity. The carrier frequency iVc, which 
is close to the atomic frequency iVa, is locked on the cavity 
mode frequency cj, while the sidebands are assumed far off 
resonance, typically atomic linewidths. By demodu¬ 

lating the light transmitted through the cavity at the FSR, 
detection of the non-linear phase response of the transmitted 
light is achieved. This phase response results in a photocur¬ 
rent that serves as a frequency discriminating error signal that 
stabilizes the laser frequency through the requirement Ai = 0. 


input intensity in the region where bistability is present 
would in principle allow the greatest degree of stabi¬ 
lization |I3l| . However, it is not practical to work in 
the bistable region since quantum and classical fluctu¬ 
ations between the semi-classical eigenmodes cause the 
system to be dynamically unstable. Therefore, one is re¬ 
stricted to working at input intensities above the bista¬ 
bility, where the achievable stabilization is orders of mag¬ 
nitude worse. Still, it was shown 13| that phase stability 
corresponding to the sub-mHz level should be achievable. 


Recently [H, an experimental effort to demonstrate 
the cavity QED system was made by probing the — 
|®Pi) intercombination line of ®®Sr atoms (i.e. 7a/27r = 
7.6 kHz) inside an optical cavity. There, however, the 
atoms were not trapped in an optical lattice, but loaded 
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into the center of the cavity using a MOT, which was then 
turned off during probing. Typical MOT temperatures 
correspond to a few millikelvin which is equivalent to a 
Doppler width of several MHz. Considering the narrow 
7.6 kHz line of the optical transition, this implies that 
motional effects will be important. 

In this paper, we extend the many-atom cavity QED 
theory of [l^ to include atomic motion, and study its 
effect on the stabilization precision. In spite of the large 
Doppler effect, the standing-wave nature of the cavity 
field induces sharp saturated absorption and dispersion 
features to appear in the considered observables. These 
features are nestled in the center of the overall Doppler 
broadened features [Mi- The stabilization that is 
achievable by utilizing these sharp features is impeded 
by multi-photon scattering processes that occur when an 
atom’s velocity matches one of its Doppleron resonances 
[m [ 2 ^ . The dependence of the stabilization on the num¬ 
ber of atoms and the temperature due to the Dopplerons 
is discussed. We demonstrate that the motion of the 
atoms causes the bistability region to disappear, so that 
no restrictions on input power are necessary to avoid the 
dynamic instability that would otherwise result. 


MODEL 


We model our system as a collection of N two-level 
atoms inside a single mode optical cavity using the quan¬ 
tum Born-Markov master equation to describe the open 
quantum system. 




+ ^ [p]; 


( 1 ) 


where. 


H = — a^ + fiTj (a^' -b a) +h'^gj{t) + (T+ a) , 

1=1 1=1 

( 2 ) 

and £ [p] denotes the Liouvillian. 

The Hamiltonian H describes the coherent evolution of 
the coupled atom cavity system in an interaction picture 
which rotates at the frequency of the cavity, and A is the 
atom-cavity detuning. The Pauli spin matrices for the 
atoms are ,tT“ and (t|, and a is the annihilation opera¬ 
tor of the cavity mode. Furthermore, p = (nPin)/{huj) 
is the classical drive amplitude, where k is the decay rate 
of the cavity, Pi„ is the input power, and uj is the fre¬ 
quency of the cavity mode. The atom-cavity coupling 
rate is gj{t) = go cos{Sjt), where go is the maximum cou¬ 
pling amplitude, and 6j = kvj is the Doppler shift in 
terms of the velocity Vj of the jth atom, and wave num¬ 
ber k of the light. 

The incoherent evolution describes the various forms 
of dissipation in this system and is described by the Li¬ 


ouvillian C [p ], 

£[p\ = —— {a^O/3 -b pa^a — 2apa^} 

N 

+ p^t - 2-^7 p^t) 

i=i 

+ ^Y{^Vp^'o-p)^ ( 3 ) 

i=i 

where p is the system’s density matrix, 7 is the spon¬ 
taneous emission rate for the atoms, and 7 p is the total 
decay rate of the atomic dipole. 

We derive Langevin equations corresponding to 
Eq. (IT|). Assuming that the classical drive rj is suffi¬ 
ciently strong, a mean-field description provides an accu¬ 
rate representation [ 2 ^. We define the mean values for 
the field a = —i (a), and for the atoms, a~ = (d”) , cr^ = 
(d^) ,(t| = (d^), which evolve according to the semi- 
classical evolution, 


N 

a = -Ka + r] + '^gj{t)a~, (4) 

1=1 

d” = - (7p + * A) a~ + gj , (5) 

= -7 + 1 ) - ‘i.gjit) (acr+ -b a aj ) . (6) 

In the moving frame of reference of the jth atom, the 
cavity field appears as a traveling wave, containing two 
frequencies shifted above and below the cavity frequency 
by the Doppler shift 5j (refer to Fig.[T|). 

It is convenient to approximate Eqns. (I4l-[6|) as a func¬ 
tion of the continuous variable 5 = kv, 

a = —Ka + r] + go J dSP{S) cos{6t)<7~, (7) 

d“ = — ( 7 p -b iA) a~ -b go cos{St)aa^, ( 8 ) 

d^ = —7 (cr^ + 1) — 2go cos{St) (acr'*' -b a*cr“) , (9) 

where P{S) is the Maxwell velocity distribution of width 
(5o, which is related to the temperature by the equiparti- 
tion theorem. 

To solve this problem that intrinsically contains a bi- 
chromatic drive, we proceed in two ways. Our first ap¬ 
proach is to numerically integrate Eqns. dllEl), partition¬ 
ing the integral into finite-size velocity bins. The veloc¬ 
ity partition must be chosen with care, since the system 
exhibits Doppleron resonances, which have a strong de¬ 
pendence on the atomic velocity. Specifically, at lower 
velocity, more resolution in the partition is required. 

Our second approach is semi-analytic, and involves a 
Floquet analysis [ 2 ^, in which we expand (T“, (T+, 





3 


and in terms of their Fourier components, 
a- 

I 

I 

= ( 10 ) 

I 

where Xi\x 2 \ and are the amplitudes of the Ith 
Fourier component. Upon substitution of Eqns. m into 
Eqns. dZHH), equations for the amplitudes are found: 


ii'' = - (*(^ + + 7p) a:i'^ + ^ + 4 > 

( 11 ) 

4^ = (*(^ + - 7p) 4^ + + 4”^^) > 

( 12 ) 

4 ^ = -7'5z,o - {US + 7 ) 

—go + ctX2~^'^ + Q;*a:4^^) • 

(13) 

In order to find a steady state solution, we set the 
time derivatives in Eqns. (Hiiini) to zero, and substitute 
Eqns. (fTTHT^ into Eq. (fT^ . yielding. 


0 = 7<^i,o + aix^3^ + + bix^3 ~^^, (14) 

where Si^o is a Kronecker delta, and 

1111 


ai = US + ^ + 




Qi+i Pi+i Qi-i Pi-i 




<3;-i Pi-i 


(15) 

(16) 


where. 


resonance corresponds to the terms of order //2 in the 
Floquet theory 2ll 22|. This correspondence is why only 
even values of I can couple into Eq. Cl¬ 
in steady state, Eq. 0 simplifies to, 


K 2k 


J dSP{S) (^x[~^'> + x[^'>y (19) 


Since x-^^^ and x\ depend on a, the self-consistent field 
amplitude a that solves Eq. CD is found numerically by 
applying Newton’s method for root finding (^ . 

We have seen excellent agreement between the two pre¬ 
viously described solution methods, and for the remain¬ 
der of the paper, focus our attention on the Floquet so¬ 
lution, which most transparently illuminates the under¬ 
lying physics. 


DISCUSSION OF STEADY STATE SOLUTIONS 


We first consider the lowest order solution to Eq. m, 
by truncating at Z = 0 , which means we have not included 
higher order Doppleron processes. We have verified that 
this solution displays the correct physics qualitatively by 
comparing to higher order solutions that are truncated 
at increasing values of 1. We define scaled intracavity 
and input field amplitudes x = af^/ng, y = r]/{Ky/no), 
where no = (77p)/(4go) i® saturation photon number. 
Eq. (fT^ then becomes. 


y=x{ 1 -b 


NCo 


dSP{S) 


1 - i{A + S)/"fp 




1 - Z(A - (5)/7 


1 + 





4 ^( 1 +? 


t± _ 7p-KA±<5)’’ 


( 20 ) 


where Co = and = -ttttwxw- 

It is interesting to consider the relation corresponding 
to Eq. (I 20 I) for a ring cavity system that has a field trav¬ 
eling in only one direction, assuming equal intracavity 
power. The corresponding expression is given by, 


Pi = i{l5 -b A) -b 7p, 

(17) 

Qi = i{lS - A) -b 7p. 

(18) 


For a given a, Eq. (HI defines a tridiagonal linear sys¬ 
tem that can be solved by truncating I at some finite 
value, and applying the Thomas algorithm for matrix in¬ 
version [ 2 ^. 

Since the atoms have motion, the condition for reso¬ 
nance between an atom and photon is achieved when the 
atomic frequency and photon frequency are offset by S. 
However, higher order multi-photon processes, known as 
Doppleron resonances, involving 2n -|- 1 photons where n 
is an integer, are also possible. The nth order Doppleron 


y = x\ 1 + 


NCo 


d6P{S) 


1 - i{A + S)/jp 


1 


2 


( 21 ) 


In the experiment described in 16|, the measured ob¬ 
servables are the cavity transmitted power T = Ix/gp 
and transmitted phase shift (j) = arg{x/y) of the intra¬ 
cavity light relative to the input light. Fig. [5] shows that 
the presence of in Eq. (1^ results in extra absorption 
and dispersive features (blue solid) around resonance in 
the transmission and phase shift, as compared to a ring 
cavity field (red dashed) where = 0 . 

These extra absorption and dispersive features are 
caused by the following: The distribution of atomic veloc¬ 
ities results in a different Doppler frequency shift of the 
























4 






10 “ 10 ‘, ,, 10 ^ 10 ^ 
lyr 


features caused by the two component saturation and 
the features caused by the single component saturation 
are clearly able to be distinguished. It can be seen in 
Fig- [21(c) that as the amount of saturation becomes large, 
the central feature becomes power broadened, but is still 
identifiable. Fig. [21(d) shows the intracavity intensity for 
a given input intensity with the values used in Figs. [H 
(a-c) labeled by black dots. 


SHOT NOISE LIMITED LASER STABILIZATION 

For laser stabilization, the central part of the phase 
response close to the atomic resonance serves as an error 
signal, see Fig. |21 (c), and allows for the generation of a 
feedback signal to the laser frequency. Through the pho¬ 
todetector, the error signal is converted to a measurable 
photocurrent. Any photocurrent measured by the photo 
detector should, in principle be zeroed by an ideal feed¬ 
back loop of infinite bandwidth to the laser frequency. 

To determine the potential phase stability that could 
be achievable using our system, we consider the shot noise 
limited stabilization linewidth assuming a strong local 
oscillator, as derived in 0, 


FIG. 2. (color online). Development of the extra absorption 
and dispersion features in the transmission (T) and phase 
shift (0), as the input intensity is increased. Blue (dark gray) 
solid curves are for a standing wave cavity, red (light gray) 
dashed are for a traveling wave cavity. For all plots, NCo = 
600 and So/jp = 260, which in the case of ®®Sr, corresponds 
to a temperature of ~ 15 mK. (a) |i/p = 5 x 10~^, (b) 
\y\'^ = 6 X 10^, (c) \y\'^ = 9 x 10^; The inset is zoomed in 
to emphasize the central features, (d) Input vs intracavity 
intensity at resonance. The black dots label the input and 
intracavity intensities of (a) (b) and (c), and the dashed line 
is |?/|^ = \x\‘^ for reference. Note that there is no bistability. 


light for each atomic velocity class, so that each velocity 
class will be resonant at a different detuning A. When 
|A| <C 7 p|a;|, the resonant velocity class of the atoms is 
interrogated by both components of the standing wave 
field, whereas when |A| 7p|a;|, the resonant velocity 

class of atoms is interrogated by only one component of 
the standing wave field. Thus, there is an increased sat¬ 
uration in atomic absorption, with a corresponding sat¬ 
urated dispersive feature, for |A| <C 7p|a;|. These sharp 
features are absent from the traveling wave cavity situa¬ 
tion where there is only one propagating field. 

As seen in Fig. [2] (a), for low input intensities, there 
is no atomic saturation at any detuning. In Fig. [2] (b), 
the atoms in the velocity class around resonance are sat¬ 
urated by both components of the field, and the velocity 
classes away from resonance are saturated by only a sin¬ 
gle component of the standing wave field. Therefore, the 


Ah' = 


hu 


Cn 


SirePf 




( 22 ) 


V / 


where Psig is the signal power, e is the photo-detector ef¬ 
ficiency and (^) is the dimensionless phase slope at res¬ 
onance A = 0. In the case that FM spectroscopy, such as 
NICEOHMS [l^ , is used for the detection of the absorp¬ 
tion or cavity transmitted phase, Eq. must be mod¬ 
ified as follows. In general, in a configuration where side¬ 
bands are applied at the free spectral range of the cavity, 
Eq. (I22|) must be multiplied with (1 + Psig/2Psideband), 
where Psideband is the sideband power. 

Eq. is the smallest, and hence the phase is most 
stable, when the product of the slope around resonance 
and the intracavity intensity is as large as possible. The 
optimal input intensity, which allows this product to be 
as large as possible, is the value used in Fig. |2] (c), and 
is labeled by the black dot at \y\‘^ 10^ in Fig.|^(d). 

To achieve a quantitative agreement between theory 
and experiment 16|, higher orders in I in Eq. (1191) must 
be included. These higher order terms correspond to 
Doppleron resonances, i.e. multi-photon scattering pro¬ 
cesses between the atoms and cavity mode. 

To study the importance of these higher order 
Dopplerons, we calculate the linewidth from Eq. (1221) at 
the optimum input power while varying the order of I at 
which the truncation occurs. 

Fig. m shows the dependence of the linewidth at the 
optimum input intensity on the order of I at which the 
truncation occurs for 3 different sets of parameters. The 
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FIG. 3. (color online). Linewidth at the optimum \y'\^, which 
gives the smallest linewidth, as a fnnction of included orders 
of I, normalized by the I = 0 linewidth. Blue circles: NCo — 
600, 5o/7p = 260, \yf = 1270. Red squares: NCo = 6000, 
<5o/7p = 260, \y\^ = 7700. Black diamonds: NCq = 600, 
5ohp = 80, \y\^ = 850. 


first set, shown in blue circles, converges by I = 12. The 
linewidth calculated with up to / = 12 included before 
truncation is around 5 times larger than the linewidth 
with only I = 0 included. This shows that Doppleron 
effects are crucial to include for a correct quantitative 
analysis of this system. 

Shown in red squares, we increase the number of atoms 
by a factor of 10, and again calculate the linewidth as a 
function of the order of I at which the truncation occurs. 
As a result of N being increased, the optimum value of 
|t/p is also increased. Convergence occurs around I = 
16. Now, the difference between the converged linewidth 
and the linewidth with only I = 0 included has increased 
by a factor of around 2. This demonstrates that as N 
is increased, higher order Dopplerons play an increasing 
role. 

We also decrease the temperature by a factor of 10, 
and again calculate the linewidth at the new optimum 
value of |j/p, as shown by the black diamonds. Even 
though the optimum |?/p occurs at a lower value, there is 
still an increase in the difference between the converged 
linewidth and the I = 0 linewidth. Convergence occurs 
around I = 12. This shows that as the temperature is de¬ 
creased, higher order Dopplerons also play an increasing 
role. 

We next study the effect of optical bistability, and its 
effects on the optimum input intensity. As can be seen in 
Fig. m (a), when So/jp = 0 there is a bistability in input 
(|t/p) vs intracavity (|a;p) intensities. As the temperature 
is increased, the bistability becomes less pronounced, un¬ 
til eventually it disappears entirely, which can be seen in 
the So/jp = 30 curve. 

The slope at resonance, which can be seen in Fig.|4](b), 
has very a different behavior for the zero and non zero 
temperature cases. For do/7p = 0, the slope is largest on 
the lower branch, decreases slightly in the bistable region, 
and then drops exponentially to zero as the intensity is 
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FIG. 4. (color online), (a) Scaled input {\x'^) vs intracavity 
(lyl^) intensities for several temperatures with NCo = 800 (b) 
Slope at resonance for several temperatures with NCo = 800 
(c) Stabilization linewidth as a function of NCo for several 
temperatures. For (5 o/7p = 0, black dashed is calculated at 
the fixed input intensity |i/P = 4ACo; black solid is calculated 
at the fixed input intensity |i/|^ = (AC'o)^/4. For So/'jp 7 ^ 0, 
the linewidth is calculated with the input intensity fixed to 
the optimum value that gives the smallest linewidth. In the 
case of ®®Sr, (5 o/7p = 150 corresponds to a temperature of ~ 
5 mK. 


increased on the upper branch. When the temperature 
is increased, the bistability disappears, and a dispersion 
feature with negative slope appears. Then the slope goes 
to a maximally negative value before power broadening 
eventually causes it to increase to zero. 

The disappearance of the bistability is important for 
achieving the highest possible degree of stabilization. In 
the zero temperature case, the optimal combination of 
\x\^ and ^ to give the smallest Aiz occurs on the far left 
side of the upper branch, when the input intensity has 
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Transision 

A 

7/27r 

F 

N 

NCo 

Temperature 

P. 

in,opt 

Au 

^So ^-^Pi 

556 nm 

182 kHz 

250 

1000 

2.5 • 10 '^ 
5.0 • 10'^ 

374 

2991 

6.5 mK 

128 fiW 

312 

3.3 Hz 

207 mHz 

'‘°Ca ^So -^®Pi 

657 nm 

400 Hz 

250 

1000 

2.5 • 10 '^ 
5.0 • 10'^ 

522 

4176 

1.7mK 

5.5 nW 

2.2 nW 

322 mHz 

20 mHz 

^So ^^Pi 

457 nm 

34 Hz 

250 

1000 

2.5 • 10 '^ 
5.0 ■ 10 '^ 

253 

2021 

3.0 mK 

0.5 nW 

0.2 nW 

8 . 1 Hz 

500 mHz 

®®Sr ^So ^®Pi 

689 nm 

7.6 kHz 

250 

1000 

2.5 • 10 '^ 
5.0 ■ 10 '^ 

574 

4593 

3.0 mK 

47 nW 

84 nW 

102 mHz 

6.8 mHz 


TABLE 1. Shot noise limited linewidths Ai/ estimated theoretically for a number of ^So —^-^Pi intercombination lines at 
experimentally realizable parameters: transition wavelength A, natural linewidth 7 , atomic sample temperature as reported in 
[27 h^ and optimal cavity input power Pin,opt- The critical temperature at which bistability disappears for NCo = 4593, the 
largest value of NCo here considered, is 60 fiK, which is well below any temperatures considered. The ratio of the carrier and 

P ■ 

sideband power is chosen to = 1 in all cases. Throughout the calculation we have assumed a prestabilized probe laser 

decoherence of yiaaer/Syr = 2.0 kHz. The shot noise limited linewidths are calculated for systems with empty cavity finesse of 
F = 250 and N = 2.5 • 10^ atoms overlapping the cavity mode, and an improved system with F = 1000 and N = 5.0 ■ 10^. 


the fixed value |yp = AN Co, labeled with the black dot 
in Fig. [4] (a) and (b). However, this value of intensity 
is in the bistable region, so the system is dynamically 
unstable when the full quantum dynamics are accounted 
for. If the tunneling rate between the different branches is 
not small, we are constrained to work at input intensities 
that are above the bistable region where the achievable 
frequency precision is much less. The grey dot in Fig. 0] 
(a) and (b) is at the far right side of the bistable region, 
at the fixed input intensity value |yp = (A^Co)^/4. The 
stabilization linewidth is orders of magnitude worse here, 
because the slope is so small. Nonetheless, working above 
the bistability allows for shot noise limited stabilization 
line widths of ~ 1 mHz [l3| . 

Fig. [4] (c) shows the stabilization linewidth as a func¬ 
tion of NCo for several temperatures. In the So/jp = 0 
limit of zero temperature, the black dashed curve corre¬ 
sponds to the input intensity corresponding to the black 
dot, and the black solid curve corresponds to the in¬ 
put intensity of the gray dot. When the temperature 
is increased sufficiently, the bistability disappears. Then, 
there are no longer any regions of dynamic instability 
and therefore no restrictions on input intensity. Each of 
the 60 /"fp ^ 0 cases in this plot are calculated at their 
respective optimal input intensity, and stop at a critical 
value of A^Cq, where the bistability reappears. It follows 
as reciprocal information that for a given value of NCo, 
there exists a critical temperature at which the bistabil¬ 
ity disappears. Therefore, by using Fig. 0] (c) to note 
the values of temperature and NCo at which the bista¬ 
bility disappears for the different curves, it is seen that 
the critical temperature increases as NCo increases. 

In general, a lower temperature will yield a smaller 
linewidth. However, because one must avoid working in 
a region of optical bistability, the optimal shot noise lim¬ 
ited linewidth for certain values of NCo and 5o/lp in 


which there is no bistability can actually be smaller than 
the 60 /jp = 0 linewidth which is constrained above the 
bistability. 


ATOMIC SYSTEMS 


In order to evaluate the performance of the cavity QED 
stabilization method presented in this work, we have 
estimated the shot noise limited linewidth using equa¬ 
tion (1221) for for a number of different atomic systems. 
For each element we assume the atomic sample laser to 
be cooled close to the Doppler limit with reference to 
temperatures obtained experimentally in [illiJil, see 
Table H These temperatures are all above the critical 
temperature where bi-stability disappears by at least an 
order of magnitude. 

A realistic atom number overlapping the cavity mode 
has been estimated to be Y = 2.5 • 10^ and has been 
used for all the elements. For several elements, larger 
atom numbers have been reported in the literature which 
makes our estimates somewhat conservative [Ulil-Ei. 
Cavity dimensions are based on the experimental values 
obtained in 16|. The cavity waist diameter is chosen 
to be 1.0 mm and the empty cavity finesse to be F = 
250 and F = 1000, corresponding to an empty cavity 
decay rate of k/2tt = 2.0 • 10® Hz and k/2'k = 0.5 • 10® Hz 
respectively. In each case, the laser decoherence has been 
assumed to be Tjaser “ ’ 2.0 kHz. For each element in 

Table E we have used the laser power that corresponds 
to the minimum shot noise limited linewidth. 

Whereas stabilizing a probe laser on the ^Sq — 
transition of ^^Mg may be slightly more experimentally 
challenging compared to the other elements in Table E 
the stabilization of a probe laser on the ^Sq —transi¬ 
tion of ®®Sr with cavity finesse F = 250 and N = 2.5 • 10^ 
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promises shot noise limited linewidth of 107 mHz. This 
linewidth is already comparable with the liiiewidths of 
the current state-of-the-art frequency stabilized lasers 

11 [i3- 

The shot noise limited linewidth may be reduced fur¬ 
ther by increasing the value ot NCo. This can be achieved 
by increasing the atom number or the cavity finesse. The 
shot noise limited linewidth, Ah', is inversely proportional 
to the square of the atom number according to equa¬ 
tion (| 22 ]) . as the phase slope at resonance, depends 
linearly on the atom number N. A similar scaling can 
be achieved by increasing the finesse while keeping the 
classical drive amplitude rj constant. For ®®Sr, the laser 
linewidth may be reduced to Ah' = 6.8 mHz by only in¬ 
creasing the cavity finesse to F = 1000 and increasing 
the atom number by a factor of two to IV = 5.0 • 10^, see 
Table H 

A cavity imposes a delay on the error signal possibly 
degrading the feedback to the laser. Generally, the op¬ 
tical cavity effectively acts as a low pass filter with a 
cut-off frequency given by fcut = c/{‘iLF) [s^, where c 
is the speed of light, L is the cavity length and F the 
finesse of the cavity. In the case studied in this paper 
we find the cavity belonging to the so-called bad cavity 
limit, where the cavity linewidth Tc is significantly larger 
than the the atomic linewidth 70 . For a typical cavity of 
length 10 cm, operated in this domain, the finesse ranges 
from 100-1000 and corresponds to cut-off frequencies in 
the MHz range. For practical purposes this will not pose 
any significant delay in the servo loop. 

The experimental proof-of-principle demonstrated in 
0 is operated in a cyclic fashion, on a time scale similar 
to that of state-of-the-art optical lattice clocks. In these 
clocks, the preparation of cold atoms takes of order 0.5 - I 
seconds, before they are interrogated by the clock laser. 
This leaves the probe (clock) laser uncorrected during 
the dead time period. All such systems suffer from the 
so-called optical Dick effect [3lJ, where the frequency 
noise of the clock laser is aliased into the sampled signal 
and ends up contaminating the clock stability. Possible 
ways to decrease the Dick effect is to increase the probe 
period by transferring the atoms to an optical lattice or 
by reducing the initial clock laser instability. With a 
laser noise level similar to that achieved by optical lattice 
clocks, and a duty cycle of about one second, we would be 
limited to stability at the 10“^® level [s^. Compared to 
the estimated shot noise limited line widths presented in 
Table U this corresponds to the low finesse cases. Push¬ 


ing EH limit, as shown by the slightly higher finesse of 

1000 , would require the atoms to be loaded into an opti¬ 
cal lattice, or require a reduction of the duty cycle by a 
fact or 10 to 100 ms in the exporimont prosontod in jl 6 l |. 


rather than atoms prepared in a MOT and then released. 
Using experimentally achieved numbers from 3^, we 
have estimated signals comparable to that achieved in 
0 , but now allowing for a continuous interrogation. 
This seems to be a very tantalizing and promising ap¬ 
proach for these systems. 


CONCLUSION 

We have seen that thermal atoms in a standing wave 
cavity field exhibit additional phenomena that were 
not observed when considering a frozen arrangement of 
atoms. Specifically, when the detuning of the laser and 
atoms is less than the power broadened linewidth, the 
system interacts with both components of the standing 
wave field. This causes new absorptive and dispersive 
features in the observables, which are the features that 
can be used as an error signal for frequency stabilization. 
Multi-photon scattering processes due to Doppleron reso¬ 
nances cause the stabilization linewidth to increase. This 
effect becomes more dominant as the collective cooper- 
ativity NCq is increased, and as the temperature is de¬ 
creased. A system with sufficient NCq and no atomic 
motion will exhibit optical bistability. Atomic motion 
may cause this bistability to disappear. When optical 
bistability occurs, one is generally restricted to using an 
input intensity that lies outside of the bistable region 
in order to prevent hopping between semiclassical solu¬ 
tions. When there is no optical bistability, there are no 
dynamically unstable regions, so that no such restrictions 
on input intensity are necessary, allowing the optimally 
smallest stabilization linewidth to be achieved. 
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